Effects of Inflation on a Cosmic String Loop Population 
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We study the evolution of simple cosmic string loop solutions in an inflationary universe. We show, 
for the particular case of circular loops, that periodic solutions do exist in a de Sitter universe, below 
a critical loop radius RcH = 1/2. On the other hand, larger loops freeze in comoving coordinates, 
and we explicitly show that they can survive more e-foldings of inflation than point-like objects. 
We discuss the implications of these findings for the survival of realistic cosmic string loops during 
inflation, and for the general characteristics of post-inflationary cosmic string networks. We also 
consider the analogous solutions for domain walls, in which case the critical radius is RcH = 2/3. 



I. INTRODUCTION 

Topological defects are unavoidably formed at cosmo- 
logical phase transitions [H, Q . Studying their physical 
properties, evolution and cosmological consequences is 
therefore mandatory for a proper understanding of the 
early universe. The last three decades have seen dramatic 
progress in this task (see [3] for a review), but significant 
knowledge gaps still exist. The aim of the present paper 
is to eliminate one of these gaps. 

Defect-forming phase transitions often occur near or 
at the end of inflation. Moreover, it is possible that var- 
ious stages of inflation occur, with defects being formed 
in between them 0]. It is therefore important to under- 
stand the effects of inflation on the defects, as well as 
to quantify their ability to survive any inflationary peri- 
ods that occur after they form. The effects of inflation 
on the internal (microscopic) structure of defects have 
been studied in which shows that those with thick- 
ness 5H > 1/^/2 get smeared by expansion, while those 
with smaller thickness survive. Effects on macroscopic 
(cosmological) scales are well known for the case of long 
string networks [1, Q , but this is not so for the loop pop- 
ulations, despite the fact that they are known to contain 
a total amount of energy which is cornparable to (if not 
greater than) that in the long strings [1, 0] • 

Here we study the effects of inflation on specific 
loop solutions, and consider both their microscopic and 
macroscopic evolution. We also discuss the consequences 
of our findings for cosmological scenarios involving string 
networks. Incidentally, we note that loops in a flat 
anisotropic universe were studied in where it was 
shown that the anisotropy of the background has an ef- 
fect on the loop motion; the reasons for this will be- 
come clearer in what follows. We will also present a brief 



•Electronic address: "ppavelin@fc.up.ptl 

^Electronic address: C. J.A.P.Martins@dai ]itp.cam.ac.uk| 

■f Electronic address: |E.P.S.Shellard@damtp.cam.ac.ukl 



analysis of the analogous solutions for domain walls. We 
shall assume that the source of inflation is a perfect fluid 
with equation of state p = wp (with w < —1/3) and the 
scale factor behaves as a c>c i2/3(i+iu). jf — _i then 
a (X exp{Ht) with H being constant. 



II. 



MICROSCOPIC LOOP EVOLUTION 



The world history of a Goto-Nambu cosmic string 1 1 
can be represented by a two-dimensional world-sheet 2 1 
= x^ia"-), with a = 0,11^ = 0,1,2,3, obeying the 
usual Goto-Nambu action. In a flat FRW universe the 
line element is given by ds'^ = a?{ri){drj^ — dx^), where 
a is the scale factor, ry is the conformal time and x are 
conformal spatial coordinates. If we identify conformal 
and world-sheet times and require that the string velocity 
be orthogonal to the string direction (i.e. x-x' = OVthen 
the string equations of motion take the form 0, 0,1^ 

x + 27i (1 -x2)x= i ^^"j ,e + 2Hi2e = (1) 



where the coordinate energy per unit length, e, is — 
x'^/(l — x^), Ti. = d/a, and dots and primes are deriva- 
tives with respect to the conformal time and space coor- 
dinates. 

Now consider the evolution of a circular cosmic string 
loop. Its trajectory can be described by x(?7, 0) = 
q{r]) (sin 6', cos 0,0) with < 6 < 2tt. Let us define the 
invariant loop radius as i? = |g|a7, proportional to the 
energy of the loop (it is also useful to define the physical 
radius, r — aq) and 7 = (1 — v^)^^/^, with v = q being 
the microscopic speed of the loop. Then the microscopic 
string equations of motion ([T]), written in terms of R, v 
and physical time t, are 



dt ^ ' dt 



fiv 



^ 2Hv) (2) 



where H = H/a is the Hubble parameter and f(v) 
'ysign{-q). 



2 



w=-1, rH=0 w=-1, rH=0.3 




-1.5 -1 -0.5 0.5 1 1.5 20 40 60 80 100 120 

r/r^ r/r. 

FIG. 1: Phase space evolution of four circular loops (with 
different initial sizes) in an Einstein-de Sitter universe. 

We start by studying the evolution of four different 
initially static circular loops (with a range of initial sizes) 
in a de Sitter universe {w = —1). In Fig. [T]we plot their 
phase space diagrams (r/ri,v): their motion is periodic 
as long as the loop is small enough compared with the 
Hubble radius, that is 

r,H <R,H=^. (3) 

This happens because if H, is a constant then Eqns. ([2]) 
are invariant with respect to the transformation v —v, 
f — /. This symmetry of the circular loop equations of 
motion implies that if the loop is able to collapse, which 
is the case for < Rc (as in this regime the effect of 
curvature dominates over the maximum Hubble damping 
effect), then the loop motion will necessarily be periodic. 
One can therefore say that the motion of a loop with 
initial radius smaller than the critical radius given by ([3]) 
never becomes dominated by the background cosmology. 
Nevertheless, smaller loops than this critical size do get 
affected by it. This can be easily seen, for example in 
Fig. [2] where we plot the loop evolution as a function 
of cosmic time: the distortion of the curves caused by 
the very rapid expansion of the universe is evident. In 
particular the loop period becomes larger as we increase 
their size, approaching infinity when — > i?c. 

For periodic solutions {d{\nR)/dt) — and conse- 
quently it follows from Eqn. ([2]) that {v"^) = 1/2 if 
is a constant. Here the brackets denote a time av- 
erage over one period. This is expected since the en- 
ergy density in sub-critical periodic loops should evolve 
as matter [p (x a"^) assuming a negligibly small loop 
production from the long string network. Since the cos- 
mic strings equation of state is w = p/p — (2(t;^) — l)/3 
then (u^) = 1/2 is just what we need for w = 0. 
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FIG. 2: The time evolution of the loop radius and velocity, 
for the loops of Fig. [1] The solid and dashed lines correspond 
to the loop size and velocity, respectively. Time is in units of 
the initial loop length. 

There is also a stationary solution of Eqns. 12]) with 
fixed physical radius and velocity 

Hr = l/V2, f = 0. (4) 

This may appear to be static but it's only stationary — 
it is a contracting loop standing still against the Hubble 
expansion. Note that in comoving coordinates one has 
Tiq — 1/^/2 and v ^ q = — 1/\/2. Above we assumed 
the loops to be initially static {vi = 0), which does not 
include this limiting case. However, with = and 
ri — Rc, our solutions asymptote to the stationary so- 
lution. For larger loops with > i?c the motion is not 
periodic and the loop will freeze in comoving coordinates. 
Specifically, the loop will asymptotically behave as 

q — const. -4=^ R oc a , v oc . (5) 

This behavior is similar to the corresponding scaling law 
for the long string network, discussed in . 

The situation is more complicated if w ^ —1 and/or for 
non-circular loops: then the loop motion will no longer 
be periodic. An example can be seen in Fig. [3l Given 
that for w > —I the Hubble radius will be a growing 
function of time, the de Sitter case will actually be the 
best possible situation for loop survival in the absence of 
phantom behavior. In all other cases, loop survival will 
be more difficuh and for R < R^ ^ O.bR-^ = 0.75(l+w)i 
the loops slowly disappear by decaying into gravitational 
radiation. Although we only study very simple loop so- 
lutions we expect these results to approximately hold for 
realistic loops chopped off by a cosmic string network. 
Most of these are produced on scales smaller than the 
horizon, typically at the scale of the string correlation 
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FIG. 3: Phase space diagram and time evolution of a circular 
loop with near-critical size in the case w — —0.5 (top) and of 
an initially elliptic loop with a major to minor axis ratio of 2 
(bottom). Solid and dashed lines respectively correspond to 
the loop size and velocity, and time is in units of the initial 
loop length. For the elliptic loop there are two sets of curves, 
corresponding to projections along the major and minor axes 
(respectively x and y). 



FIG. 4: Comparing the full (numerical, Eqns. (O) and ap- 
proximate (analytic, see Eqn. ([9])) VOS model solutions for 
the evolution of large loops during an inflationary epoch. 
These are respectively shown by solid and dashed lines, for 
each set of initial conditions. Each panel corresponds to a dif- 
ferent value of the equation of state parameter w, and within 
each one the four sets of lines correspond to initial loop sizes 
e — 10, 5, 3, 1, respectively from top to bottom. 



length if there is no significant amount of small-scale 
structure (which we would expect after a few e- folds), 
or even on much smaller scales if small-scale wiggles are 
present. However such very small loops will be irrele- 
vant as far as the ability of cosmic strings to survive an 
inflationary era is concerned. 



III. AVERAGED STRING EVOLUTION 

The averaged properties of a cosmic string network can 
be accurately described in the context of the velocity- 
dependent one- scale (VOS) model It describes 
the string dynamics in terms of two macroscopic param- 
eters, the string root-mean-squared (RMS) velocity, v, 
and the string correlation length, L, which is assumed to 
be the same as the string curvature radius and is related 
to the string energy density via p = n/L^ (with /i being 
the string mass per unit length). We note that we will 
be using the variable v to denote both the microscopic 
and the string RMS velocity depending respectively on 
whether the microscopic or the averaged cosmic string 
evolution is considered. 

In theis context the evolution equations for these quan- 
tities have the following form [ 



dL , 9, dv , / k(v) 



2Hv 



where c is the loop chopping efficiency and k{v) is a func- 
tion of velocity known as the momentum parameter 
One can also use the model to describe the evolution of a 
given cosmic string loop; circular loops are characterized 
by their radius R (now an averaged quantity). In this 
case the evolution equations are 



dR , dv , 

-.(l-2.^)i/E,-.(l-.^) 



k(v) 



2Hv 



(7) 



(6) 



Note the obvious similarities with ([2]). 

We can use this model to study the evolution of the 
loops larger than i?c, and compare it with the evolution of 
the long-string correlation length. We solve numerically 
the Eqns. ([7]) or ([S]) for given values of the equation 
of state w and the initial loop size or string correlation 
length, which we parametrize as RiHi = e and LiHi = e 
respectively (note that for an initially static loop Ri = 
Ti). We expect large loops to be conformally stretched, 
so the relevant quantity to plot is Cioop = Rui/Ria for 
loops of radius R, or Qong = Lai/LiU for the long-string 
correlation length, as a function of = a/a^. Examples 
can be found in Fig. 2] for loops, and in Fig. [S] for the 
long strings; the solid lines correspond to the numerical 
solutions. 

The difference between the two results is more appar- 
ent than real and has to do with the fact that we did 
not take into account an energy loss mechanism for large 
loops, and moreover R and L are not directly related to 
each other. Recall that R is directly proportional to the 
loop's energy density, while L is inversely proportional to 
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(the square root of) the energy density of long strings in 
a given volume, = ^V/ E. Therefore, the dynamical 
effect of expansion will be reflected differently in R and 
L. We can see this easily by defining an effective long 
string energy lengthscale that is proportional to the en- 
ergy, E — ^jlV/L^ = ^Leff- Substituting this definition 
in Eqn. ^ we find 



dL 



dt 



(1 - 2v^)HL 



cvN, 



(8) 



which is similar to Eqn. ([7]) except for the loop produc- 
tion term (which is now negative since loop production 
removes energy from the long strings). We have defined 
Neff = Leff/L which is the number of segments of size 
L in the chosen volume. The energy density of both large 
loops and long strings will eventually scale as instead 
of a^"^ (the case of non-relativistic point- like particles). 
Consequently both large cosmic loops and long strings 
have the ability to survive more e-foldings of inflation 
than point-like objects. 

Since we know that the string motion will be damped, 
we can use the VOS model prediction for the relation 
between the velocity and the correlation length in this 
regime, namely [Bl, [y| w ^ ko/2Hr, where fco = 2A/2/7r is 
the value of the momentum parameter k(v) in the non- 
relativistic limit V — > 0. This can then be substituted in 
the evolution equations, and the following solutions are 
obtained 



^2 

^loop 



= 1 - 



long 



kojko + c) 
2e2(l + 3w) 



gl+3u; 



1). 



(9) 



(10) 



In fact these solutions are also valid for the de Sitter 
case, w = —1. The equation of state factor is negative, 
and hence one always has C > 1 for long strings, and 
C < 1 for loops. Figs. S] and [5] also show (in dashed 
lines) this solution for comparison with the numerical so- 
lution of the full VOS model. This analytic approxima- 
tion works very well for very large loops, and it becomes 
less accurate for relatively smaller ones. These loops do 
not have non-relativistic velocities throughout their evo- 
lution: they can be mildly relativistic for the first few 
e-foldings, before they suffer enough stretching to slow 
them down, in which case the above ansatz is a good ap- 
proximation in the early stages of the evolution. On the 
other hand, the the loop solution ^ also breaks down 
when the size is very small. 



IV. INTERLUDE: DOMAIN WALLS 

We will now provide a brief discussion of the case of 
domain wall networks. At the conceptual level the anal- 
ysis is of course entirely analogous to the case of cosmic 
strings, but the different co-dimension of the defects will 
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FIG. 5: Comparing the full (numerical, see Eqns. ^) and 
approximate (analytic, see Eqn. ([10}) VOS model solutions 
for the evolution of long strings during an inflationary epoch. 
These are respectively shown by solid and dashed lines, for 
each set of initial conditions. Each panel corresponds to a dif- 
ferent value of the equation of state parameter w, and within 
each one the three sets of lines correspond to an initial long- 
string correlation length given by e = 0.5,1,5, respectively 
from top to bottom. 



introduce some quantitative differences, and determin- 
ing where and why they emerge is helpful in the overall 
understanding of the problem. 

Thin spherical domain walls in Minkowski space have 
a conserved wall invariant area given by Sw = 4,Trjr'^ 
where as before 7 = (1 — t)^)"^/^, v is the domain 
wall velocity and r is the coordinate radius of the do- 
main wall. This implies that — j^/^r is a con- 
served wall invariant radius — the invariant circular cos- 
mic string loop radius is given by Rs — ^r. Now 
consider the evolution of a spherical domain wall in a 
flat FRW universe. Its trajectory can be described by 
x{ri,6) ~ q{rj) (cos cos 0, cos 6* sin 0, sin 0) with < (j) < 
27r, Q < 6 < TT. Then the microscopic wall equation of 
motion can be written as 



^ = (1--. )HR^^-^{l-v) 



9{v) 

Ryj 



3Hv 



(11) 

where now g{v) = 2^^^'^sign{—q). The differences be- 
tween the Hubble damping terms in the string and wall 
equations are related to the fact that the domain wall 
momentum per unit comoving area is proportional to 
in the case of an infinite planar domain wall (so that 
fi„7 oc a~^) while in the case of an infinite straight string 
it is the momentum per unit comoving length that is pro- 
portional to a^^ (so that vj oc a^^). 

As for circular loops, the motion of spherical domain 
walls in de Sitter space is periodic as long as the domain 
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wall is small compared with the Hubble radius, that is 

nH < RcH = ^ . (12) 

Again, notice that if iJ is a constant then Eqns. (fTTj) are 
invariant with respect to the the transformation v —v, 
g —g which implies that if the wall is able to col- 
lapse then the wall motion will necessarily be periodic. 
For periodic solutions {d{hi R) /dt) = and consequently 
(v^) = 2/3 if is a constant. This happens because 
the energy density in sub-critical periodic domain walls 
should evolve as matter (with p cx a^'^) assuming a neg- 
ligible energy transfer from super to sub-critical domain 
walls. The equation of state parameter for a domain 
wall gas is w = p/p = —2/3 + (w^) and consequently 
(v^) = 2/3 yields w 0. 

One can similarly find a stationary solution of equa- 
tions of motion with fixed physical radius and veloc- 
ity, representing a contracting domain wall standing still 
against the Hubble expansion. It has the form 

Hr = ^/3/8, f = (13) 

where r = aq is the physical radius, and in comoving co- 
ordinates Tiq = v^S/S, V = q = —^2/3. These are the 
analogues of Eqn. For walls with > Rc the mo- 

tion is no longer periodic and the domain wall will soon 
become frozen in comoving coordinates. Specifically, the 
domain wall will asymptotically behave as 

q = const. R oz a , v cc , (14) 

and again if H is not constant and/or for of non-spherical 
domain walls the wall motion will not be periodic. The de 
Sitter case is again the best possible context for domain 
wall survival (assuming w > —1). In all other cases, 
domain wall survival will be more difficult and we expect 
that domain walls with R < Rc will decay into radiation. 
Finally, notice that the energy density of large domain 
walls will eventually scale as instead of (the case 
of non-relativistic point-like particles) or (the case 
of long strings) . Consequently cosmic domain walls have 
the ability to survive more e-foldings of inflation than 
point-like objects or strings. 

V. CONCLUSIONS 

We studied the evolution of simple cosmic string loop 
solutions during an inflationary era, showing that large 



loops and long strings can survive more e-foldings of infla- 
tion than point-like objects. In a de Sitter universe, cir- 
cular loops below an initial invariant radius Rc — H^^ /2 
(assuming Vi = 0) have periodic motion, while larger 
loops freeze in comoving coordinates. Realistic small 
loops will simply decay into gravitational radiation while 
large loops are stretched by the inflationary expansion. 
We have also shown that loop solutions with = Rc will 
asymptote a stationary solution with constant velocity 
and invariant loop radius. These are very special loop 
solutions which we expect to be uncommon even in the 
context of string evolution in an inflationary era where 
the small scale structure can be significantly reduced. 
For non circular loops or varying H the loop motion is 
no longer periodic even in the case of small loops. 

Our results should provide a very good approximation 
to the overall behavior of realistic loops produced by a 
string network, and are also relevant for the final distri- 
bution of cosmic strings loops nucleated during an infla- 
tionary era which usually assumes that loops are simply 
stretched by the expansion. It has been shown ^ ^ that 
a cosmic string network can survive up to 60 e-foldings of 
inflation: even though it is pushed outside the horizon, 
its evolution once inflation ends ensures that it will be 
back inside the horizon in time to have observable conse- 
quences. Large loops produced during the final stages of 
inflation must therefore be included in any quantitatively 
realistic study of these scenarios. A detailed analysis of 
these issues is left for future work. 
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